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Abstract 

We study the asymptotic behaviour of isotropic and homogeneous universes in general 
scalar-tensor gravity theories containing a p — ~p vacuum fluid stress and other sub- 
dominant matter stresses. It is shown that in order for there to be approach to a de 
Sitter spacetime at large 4-volumes the coupling function, a;(0), which defines the scalar- 
tensor theory, must diverge faster than |(/)oo — for all e > as </> ^ 0oo 7^ 
for large values of the time. Thus, for a given theory, specified by Lij{(j)), there must 
exist some 4>oo G (0, 00) such that w — s- cx) and uj' juP'^'^ ^ as — > </)oo in order for 
cosmological solutions of the theory to approach de Sitter expansion at late times. We 
also classify the possible asymptotic time variations of the gravitation 'constant' G(t) 
at late times in scalar-tensor theories. We show that (unlike in general relativity) the 
problem of a profusion of "Boltzmann brains" at late cosmological times can be avoided 
in scalar-tensor theories, including Brans-Dicke theory, in which (/> — > 00 and ^ o(0^/^) 
at asymptotically late times. 

PACS: 98.80.Jk, 04.50.Kd, 95.30Sf, 98.80-k 

1 Introduction 



The cosmological approach to de Sitter space-time at large expansion times has assumed a 
double importance in the study of the universe. It provides a description of the evolution 
of the early universe after the prolonged gravitational influence of a slowly-evolving inflation 
field, whose energy density remains approximately constant for a significant interval of time. 
It also provides a good description of the late-time expansion dynamics of the universe in the 
presence of a dominant source of dark energy. We know that an exact de Sitter solution arises 
when the effective equation of state has the form p + p = 0, where p is the total isotropic 
pressure, and p is the total matter density. This can arise because of the presence of a single 
vacuum 'fiuid' with this equation of state, as first suggested by Lemaitre in 1933 [1], and 
employed in the construction of the steady state universe by Hoyle [2] and McCrea [3], and 
in the deduction of its stable asymptotic behaviour by Hoyle and Narlikar [4]. It could also 
result from the presence of an imperfect fluid which possess an effective equation of state 
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of this form, as is the situation with bulk viscosity [5] , or from the presence of higher-order 
curvature terms in the gravitational Lagrangian beyond those sufficient to generate general 
relativity [6]. The correspondence between the situation in general relativity and such higher- 
order gravity theories can be understood in terms of the conformal equivalence of the two 
theories [7] . The situation becomes more complicated when other theories of gravity [8] , which 
generalize Einstein's theory, are considered because the source of the field equations changes 
and a p+p = fluid stress no longer results in a simple de Sitter space-time. The most familiar 
example of this sort is the zero-curvature Brans-Dicke universe [9], where a p+p = fluid is no 
longer equivalent to the presence of an explicit cosmological constant term in the gravitational 
Lagrangian, and as such does not produce an asymptotic approach to the de Sitter metric 
[10, 11, 12]. Instead, power-law expansion of the isotropic expansion scale factor results, 
with the power becoming infinitely large as the theory approaches general relativity. This 
difference arises because, even when p is constant, the Gp term in the Priedmann equation 
falls in time, as t~'^ , due to the fact that the gravitational 'constant' varies as G oc t~'^ 
[13]. Higher-order gravity theories with Lagrangian contributions that are 0{1/R) also allow 
late-time acceleration to occur in a variety of different ways in the asymptotic limit of low 
4-curvature, 

In this paper, we will first consider the approach to the de Sitter metric in a generalized 
scalar-tensor theory [14], in order to gain some understanding of what type of accelerated 
expansion is possible in the situation where the simplest form of vacuum stress, with an 
equation of state p+p = 0, dominates the expansion dynamics. In most scalar-tensor theories, 
the introduction of a vacuum stress will not result in a de Sitter spacetime. We shall be 
specifically interested in those theories which omit late times solutions in which the space-time 
approaches the de Sitter solution. We will classify such scalar-tensor theories of gravity by 
the speed of asymptotic approach to this de Sitter limit. In all such theories, G asymptotes 
towards a constant and non-zero value. In the second part of the paper, we consider the 
asymptotic evolution of scalar-tensor theories in which the gravitation 'constant' vanishes at 
late times. A general method for finding solutions of scalar-tensor cosmological models was 
given in by Barrow and Mimoso [15, 17, 18]. It is possible to adapt their generating-function 
method to study asymptotics, as was done for another problem by Deruelle et al [16], but 
we shall adopt a more direct approach. We classify the variation of G{t) into three classes 
and then use our results to show how a solution of the 'Boltzmann brain' problem can be 
obtained in a wide range of scalar-tensor cosmologies, which includes the Brans-Dicke theory 
as a particular case. 

In Section 2, we specify the general field equations for scalar-tensor theories defined by an 
arbitrary coupling function ^(0) of the scalar field. In Section 3, we classify the possible rates 
of approach to de Sitter spacetime for different behaviours of a;(0) and provide the conditions 
under which a generalized scalar tensor theory will omit a de Sitter limit. In Section 4 we 
consider the asymptotic forms for the possible evolutions of G{t) in scalar-tensor cosmologies 
when a vacuum stress is present, using the characteristic G oc behaviour of Brans-Dicke 
theory as a benchmark. We show how theories in which ^ — oo asymptotically, and the 
coupling function of the theory, io{4>)^ asymptotes to a sufficiently large value, or — > oo more 
slowly than (j)^/"^, will not suffer from the 'Boltzmann brain' problem that besets late-time 
cosmological evolution in general relativistic cosmologies. In Section 5, we summarize and 
discuss our results. 
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2 Field Equations 



Consider a generalized scalar-tensor theory of gravity described by the action: 



1 

16^ 



d x^J—g I — u){(j))- 



+ Sm{gab,i'm), 



where we have set c = h = \; Sm is the matter action, labels the matter fields, gab is the 
metric, ^ is a scalar field and co{4>) is an arbitrary function that must be specified to fix the 
gravity theory. We define the energy- momentum tensor of matter to be Tab- By varying the 
action with respect to the metric, and then with respect to (f), we obtain the field equations 
and conservation equations: 



Gab = fTab + ^ {daCl>db<P - yab{dct>f) + ^ (VaVft^ - n<P) , 

T\,a = 0. 



SttT 
3+2u>' 



(1) 

(2) 
(3) 



The familiar Brans-Dickc theory arises when the arbitrary coupling function uj{(j)) is taken to 
be a constant. Under the assumptions of homogeneity and isotropy and spatial flatness the 
metric becomes: 

dr2 



1 — kr'^ 

Here k determines the curvature of the t = const slices. We assume that the universe is filled 
with some vacuum 'dark energy' with density po and equation of state po = —po, and that 
the rest of the matter has energy density pi = poK{y) and pressure pi, where y = Ina and 
limy_^oo K{y) = limj^^oo Ky{y) = 0. The conservation equations require that po = const and 
that pi/po = — {K'/S + K). The gravitational field equations then reduce to: 



3d 



;.2 



2a 
a 



Sirpo 



(1 + Ky/3 + K) 



_3a^ 

a (f) 



3A; 



2d. 
a I 



(4) 
(5) 



3 The De Sitter Limit 

Proposition 3.1. In order for there to he an approach to de Sitter at large 4-volumes, where 
~^ the coupling function oj{4>) must diverge faster than |0oo — '^^l"^"''^ for all e > as 

4>^ (t^ooT^ 0; 4>oo < oo. 

We define a/a = H = H^/{1-F{y)), (j)^ = 87rpo/3i?^, k = jH^^, and (p = (0oo -<?!') /^oo- 
The essential Einstein equations can then be rearranged to read: 



2R, 



1 _ _ 



6(1 - ip) 



"Pyy "Py _ ^) 

+{l-Ff (K^ + 27(l-(p)e-22/), 

-{l-Ff [K-^{l-^)e-^y), 



(6) 



(7) 
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These equations are equivalent to the pair: 



2Fy_ 

1 - F 



+ 6(2F - 



+ bifiy + Qif (8) 

+(1 - Ff [Ky + 6K- 47(1 - (p)e-^y) 
(12 + 12K + 3Ky) (1 - Ff 



3 + 2a; 



(9) 



For there to be a de Sitter limit we need linij^^oo H 



linij/^oo F 



and also liniy^oo Fy 



Hoc = const which implies that 
0. It is clear from Eq. (8) that F ^ 0, Fy implies 
that (fyy + 5ipy + Qif ^ because lim^y^oo K{y) = limj^^oo Ky{y) = 0, and so ipyy, (py ,<p ^ 0. 
In order to prove the proposition 3.1, we solve Eqs. (6-9) asymptotically in the limit y — > oo for 
general F{y). We then find the asymptotic form that io{(p) must take for de Sitter spacetime 
to emerge in this limit. 
We define: 



M(y) 



1 



(1 - F)2 



76 



-2y 



)■ 



We then have: 



2F 



^ + 6(2F - F^) - (1 - Ff {Ky + 6K- ije'^y) = (1 - Ff{Ny + 6N), 



1-F 

and so Eqs. (6) and (8) become: 

(l-FfMyiy) 



(fyy 



2F, 



27(1 - Ffe-'^y 



H-FYiMy+QM) = <Pyy+{5 + 



F, 



1-F 



"Py 



+ { 6 + 47(1 - Ffe-'^y + 



(10) 



(11) 



2Fy_ 

1 - F 



if. 



We consider three cases: limj^^oo My/M = 0, limj^_,oo My/M = —q for some < g < 00 
and \iuiy^ooMy/A4 —^0. If = 7 = then these cases represent, respectively: slower than 
power-law approach (of H — H^) to de Sitter, power-law approach to de Sitter, and faster 
than power-law approach to de Sitter. 



3.1 Slow approach: My/M 

We begin by considering those cases in which My/M — >-Oasy— >oo. Ifir = 7 = this 



would imply Fy/F 

\H — Hna\ — 



and so I if — H„ 



more slowly than a ^ for any r > i.e. 



00. Asymptotically, Eq. (11) reduces to; 

6M(y) ~ (ppp + 5(pp + 6(p, 
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and so 

^ ~ M{y). (12) 

Asymptotically, Eq. (10) then gives: 

-u^ipy ~ 2 - 2^. (13) 

Now if ify/ipFy A = const 7^ as y — > cx) then ln<^ ~ + const, and we can clearly see 
that we cannot have both F — ^ and <^ — ^ as ?/ — ^ as we have required for any finite A. 
Asymptotically, (9) gives: 

—io^(f'^ — 12 ~ 6L0(py. (14) 

If limj^^o l/('^¥'2/) = 0) then 

and so a; ~ Ae~^y for some A, and hence 

nm = lim > ±00 

y^oo Uify 2/— »oo ip ipy 

by lim^^oo V'l//'/? = 0. It follows then that we cannot have limj,_>oo l/(w^j/) = and hence, 
by Eq. (13), we cannot have limy^oo ^yFy/^ = 0- 

We must therefore have limj^^oo 'P/ifyFy) = 0, and so by both Eqs. (13) & (14): 

LO -. (15) 

fy 

It follows from limj^^co ipy/p = lim^^oo My/M = that lim^^oo |<y9w| — > 00. In order to avoid 
ghosts we must have u > —3/2 and so, at late times, we must have ipy <Q so that uj — +00. 
It follows that (p — 0+. In this case the gravitational constant, G = l/(f), behaves as: 

G(t) ~Goo(l + (p), (16) 

where Goo = It follows that G{t) > Goo at finite time and that the limiting value of G 
is approached from above. 

3.2 Power-law approach: liniy^^ao My/M — —q, < q < 00 

We now consider those cases where My/M —q < as y ^ 00. This implies that M{y) = 

f{y)e~'^y where fy/ f ^ as y — > 00. If = 7 = then this equates to \H — — >■ for 
all r > g 7^ 00 and \H — H^^l — > for all r < g / 0. In these cases, Eq. (11) gives: 

((6 - q)f{y) + fy) e~^^ ~ <Pyy + 5(py + 6(p, 

which, if g 7^ 6, has the solution: 

If (co = ci = 0), (0 < g < 2) and/or (cq = and 2 < g < 3) then from Eq. (10) we have 

-qf{y)e-''' -{q + q')v = -Q^u;^'- (17) 
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12 



a; ~ 



q{6 - q)ip' 

-22/ 



(18) 



ip ~ CqC 



~ cie 



-32/ 



(19) 



(20) 



It follows that: 

If Co 7^ and q>2 then: 
and from Eq. (10) we have: 

If Co = 0, ci / and g > 3 then: 
and Eq. (10) gives: 

If g = 2 then: 

if ~ 4yM(y), 
and Eq. (10) gives Eq. (19). If g = 3 then: 

(f -3yM{y), 

and once again uj is described by Eq. (20). 

If g / 6, we have found that the behaviour of uj is given either by Eq. (18), Eq. (19) or 
Eq. (20). In all three sub-cases ipu — const as y ^ oo, which is consistent with proposition 
3.1. Additionally, we see that if we define lim.r-,00 <fy / <f = —r, then if < r < 6 we have 
ip ^ and so as t ^ 00, G{t) tends to its limiting value from above. If r > 6 then (p ^ 0~ 
and G{t) tends to its limiting value from below. 

We must still deal with the q = 6 case. If g = 6 then by Eq. (11): 

where fy/f — 0. Now we therefore have if ~ coe~^^ + cie~^^ + g{y)e~^y where g solves: 

fy ~ 52/2/ - '^9y + 125- (21) 

such that \1 fy = Q then 5 = 0. If co / or (co = and ci ^ 0) then we revert to one of the 
cases considered above and u is given, respectively, by Eq. (19) or Eq. (20). We therefore 
take Co = ci = 0. It follows then from the requirement that lim_>oo fy/f = that 



lim (gyy/f) = lim (gy/f) = lim (g/f) = 0. 

2/— >oo y—too y—^oo 



Thus, we have 
Eq. (10) therefore gives: 



lim (fyy /My = lim (fy/My = llUl (f/My = 0. 



UJ 



6/(y)e" 



So. if gy/g 0, we have 



6 

/ 
9 



fiy)giy) 



(22) 



oo, 



6 



and so oj diverges faster than l/^p. If gy/g — — s, with s < 0, we cannot have fy/ f — as 
required, and so s > 0. Thus, 



\ip(jj\ 



(6 + sf 



oo, 



and once again a; diverges faster than l/tp. Finally, if gy/g — — oo, then we write g = e ^^^^ 
where b, by oo. Now Unij/^oo {yby)~^ = - limj,_^oo ^yy/^l = 0) and so \\m.y^^ gyyg / g^ = 1. 
By Eq. (21) we also have gyy ~ fy and so: 



Qfg 



9l 



6/ 6/ 

~— ^ oo, asy— ^oo. 

Gyy Jy 



It follows that if q = 6 then lo diverges faster than 

Wc have shown that if My/M ^ —q then uj either diverges as 1/ip or, if q = 6, it is 
possible for lo to diverge faster than 1/ip. In all cases then, we have, as proposed 3.1, that u 
diverges faster than l^'l"''^'*"^ for all e > 0. 

3.3 Fast approach: My/M -oo 

We now consider those cases where My/M — — oo as y — oo. It follows from Eq. (11) that 
(p ~ coe~^^ + cie~^y + (p where 

My ~ <Pyy ^<Py^M. (23) 

If Co 7^ 0, then Eq. (10) tells us that the behaviour of u> is described by Eq. (19). If cq = 
and ci 7^ then the behaviour of lo is given by Eq. (20). In both cases to diverges as ip"^. 
We now consider the case cq = ci = so <^ = From Eq. (9), we have 

-LO^ipl ~ 12 + 2L0(Pyy. (24) 

We define J hy Jy = M and J ^ as y ^ oo, and define J = e~^^y\ It follows from 
Jy/ J — ^ — oo and J — > that g, gy oo. Now lim^^oo {ygy)~^ = ~ limj^^oo 9yy/gy = and 
so: 

= 1 + ^ 1. (25) 

Jy gy 

Since <py ^ M we have ip ^ J and so, using the above relation, we have 

{2uj + uj^ip) f. (26) 

^y 

We now defined eo by —>■ for all e < eo and —>■ oo for all e > cq. We then have: 

12^7 



LO 



(1 + eo)^^' 



Since (py/(p — oc, the above expression gives us lim.y^oo = 0. From the definition of eo 
then, we must have eo > 0. Now clearly, for all e > 0, we have 

lim - — = 0, 

J/— >oo y 
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but then 



lim — — = — lim f ^ = 0. 



y^oo y 2 y—>-oo (p^ ^ 

Thus, for ah e > as ?/ ^ oo, we have: 



Therefore, by the definition of eo we must have eo < 0, but since we also found that eo > 0, 
and it then fohows that eo = 0. 

In this fast-approach case then, we have: 

u ~ -1^. (27) 

^y 

Since v'y/v? — oo, we have that uj diverges faster than in that limit, but since we also 
found that eo = 0, we have that [(^^""^cjI — > 00 for all e > 0, and so to diverges faster than 
\ip\-^+' for all e > 0. 

We also note that for uj +00, and hence to avoid ghosts, we must have 1^ ^ 0~ as 
y 00. The effective gravitational 'constant' G(t) ^ Goo(l + therefore approaches its 
limiting value from below in the late-time limit. 

We have shown that if My/M — —00 then, as per proposition 3.1, u diverges faster than 
for all e > 0. 

3.4 Summary 

We have shown that proposition 3.1 holds in each of the three possible cases. We have 
therefore proved proposition 3.1. For there to be a de Sitter asymptote we therefore require 
that 00 diverges faster than l^oo — 01"^"*"^ for all e > as y — > 00, ^ ^ <jf)oo 7^ 0. Note that 
this is a more stringent condition than one might expect from considering linear scalar-tensor 
corrections to the metric around a spherically symmetric body. In the PPN formalism, scalar- 
tensor corrections to general relativity vanish as a; — > 00 and lo' jio^ — > 0. The latter condition 
might lead one to expect a de Sitter cosmological limit if uj diverges faster than 1 / i/|i?5>oo — 
however we have found that this is not the case. In summary: for a given theory of gravity, 
specified by ((/>), a de-Sitter limit for late-time Friedmann cosmology requires that there 
exists some 0oo £ (0, 00) such that as ^ — > ^00 

for all e > 0. If we define — r = limj^_>oo H^yl'^ then — >■ for r < 6 and — > for 
r > 6. 



4 Decaying Gravity 

In theories with asymptotically de Sitter behaviour, G{t) const 7^ in that limit. Another 
interesting limit to consider is that in which — > at late times. In particular, when the 
equation of state is p = — an exact zero-curvature (fc = 0) Friedmann solution of Brans-Dicke 
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theory with this property is known. In this solution, u) = const, a cx t'^+^Z^ and G{t) oc 
If uj > 40,000 today, as imphed by tracking data for the Cassini spacecraft [19], then the 
possibihty that the evolution of our universe is described by such a solution is not ruled out 
by observational or experimental evidence. Although the G{t) evolution is strong, and appears 
to be CO independent, it is better to examine the G{a) oc a~^/(^'^+^) evolution, where the uj 
dependence appears explicitly and the lo ^ oo general relativity limit is evident. In order to 
classify the late-time behaviour of all such solutions, we define G = Gox{y)j H = HoF{y), 
k/a^ = jHQe'^y where SttGoPq/SHq = 1; and as above, we have taken y = In a. With these 
definitions, the Einstein equations are now equivalent to: 

(i-lf-t) - id + i^)-^^-'. PS) 



-(1 + 2)^ + ^ + 2^-^- (29) 



F V 2x; V2 ; X X 

We can combine Eqs. (28) and (29) to give: 

'^Py(, Xy\.. _ Qx{l + KJQ + K) Xyy ..Xy ^Te"^^ 



2^2/ (-, Xy\ _ xKy ^ 2e 2J/ X, 



1-^ = ^ + ^ + ^-^-('- + 2)ei- (31) 



F \ 2xJ F^ F2 XX X'' 

We now divide the possible rates of change of G{t) towards its asymptotic value into three 
classes and study each in turn. 

4.1 Slow Approach: Gy/G 

If G ^ more slowly than for any r > 0, then Xy/x — at late times, and Eq. (30) 
gives 

X^F^ + FyF/3, (32) 
and so using hm^^^o Xy/x = we have: 

X ~ (33) 

We can now use Eq. (31) to find a;(0): 

F ^^_KyF\_2xr_Kyx\_ 



Fy V ^Fy j Xy V 2x2, J 

Now, assume that Xylx^y ~^ for some a as y ^ oo. This would give x ~ Sexp(ai^) for 
some B ^ 0. Given that X ^ 0"*" as y ^ oo and x ^ 0, we cannot have a > — oo. Thus we 
must have Kyx/Xy ^ as y ^ oo, and so whatever K is: 

iv ~ (34) 
Xy 

We note that as y ^ oo, a; ^ oo, and also that 

,3 
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and so 



lim — 5- = 0. 



For these solutions: 

2 GqH SH 



G{t) ~ GqF^ 



Hi Sttp ' 

Furthermore, it is clear that any scalar-tensor theory where lo ^ 00 and uj' jus^ ^ as ^ 00 
will have such a late-time solution - since once is known, the formula oo ~ —2x/Xy can 
be used to find Xy hence H{y). 

All bounds from local tests of gravity are satisfied in the limit cj — 00 and u' joj^ — 0. 
Theories with stable late-time solutions of this class are therefore physically viable. We note 
that at late times, 

however the spacetime is not not de Sitter since gravitation 'constant' is decaying i.e. G ^ Q. 
4.2 Power-law approach: MvHy^aoGy/ G = —q 

We now consider solutions where a^G — > for r < q and —>■ 00 for r > q. Eq. (30) gives: 

+ 3)((Z + 2) ~ - -1^^ ^(1 + q/2). (35) 

Since we cannot have < 0, it follows that we must have limj^^oo Py/F = ~f for some r > 0, 
and 7e~^^/x — > as y — > 0, which is certainly the case if < g < 2 and/or 7 = 0, then we 
must have oc x in the y ^ 00 limit. Therefore, in this sub-case 2r = q and 

(g + 6)(g + 2)F^ 
X ^ ■ (36) 

From Eq. (28), then we have limy—^oo uj = ujoo and 

^00 = --\> -1/2. (37) 
q I 

If |7e~^^/x| — 00, which certainly requires g > 2 and 7 7^ 0, Eq. (30) gives: 



(g + 3-r)(g + 2) 



Since g > 2, we must therefore have r = 1 and also ~ —476 ^^/(g + 2)^. For such solutions 
to exist we must therefore have 7 < 0. From Eq. (28) we have: 

uq^ jq + 2f 
^ + ^--6 

Thus, we have limy—^Qo a; = Woo where 

3 
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Finally, if g = 2 and xe^^ qq = const 7^ 0, then Eq. (30) gives r = 1 and 

Hence, we must have 3ao/2 > 7. We define 7 = ftooo, and then from Eq. (28) we have 
limj^^oo uj = LOoo, where 

_ 3(1 + 260) 3 
"'"^ - 2(3 - 26o) > "2- (^^^ 

Observationally acceptable solutions may require that uo > 40, 000 today. For Wqo > 40000 
we would need q < 5 x 10~^. 



4.3 Fast approach: \Gy/G\ 00 

The final class of solutions that we consider have G — ^ more quickly than a"*", for any r > 0. 
This implies that |x?//xl ~^ For such solutions we write x = e^^^^^ and wc require that 
g, Qy 00. This further implies that Xyy ~ Xy/X- We also write F = exp(— /(j/)). If 7 = 
then Eq. (30) gives 

{fy-gy)gy^-QeMV -g)- (40) 

We define J, hy Jy = F and J — as y — 00. We note that J ~ F^/Fy as y — 00 since 
FyyF/Fy — 1. We now make the ansatz: 



With this definition of A: 

Since 51^ = 2Jy/J and Fj, = Jyy, we have 

^ _ JyyJ \ 
9y 2J2 ~ 2- 

Eq. (40) then gives 



gl ^ 12e^f-3. (41) 
It follows that A = \/3 is required. In this case then: 

Using Eq. (28) we find: 
and so as y — > 00, 

3 , , 

w ^ Woo = (42) 

If 7 7^ 0, Eq. (30) gives f ^ y + h where: 

{hy - gy) gy ~ 47 exp(2/i) (43) 



11 



We define Ly = and so as y — oo, L — oo. Solutions to this equation require g ~ /xL. 
It follows from this that hy/gy 0, and so /i = ^J^^^. Since we must have g — oo, such 
solutions only exist is 7 < 0. In this case then: 



= ~ e-^^. (44) 



Prom Eq. (28), we find: 

and so as y 00: a; — > uj^ = —3/2. 



._ j = n /4, 



4.4 Suppressing "Boltzmann Brains" 

We have shown that solutions in which the gravitational constant decays at late times i.e. 
G — 0, require that, as G"^ ~ ^ — 00 either a; — 00 or a; — ^ lOoo = const. Generally speaking, 
the smaller limj^^oo is, the faster G decays at late times. The fastest decay occurs when 
Woo = —3/2. Observational data may require thatu; > 40,000, today [19], and so solutions 
with small values of Woo are unlikely to exist in physically viable theories. Theories with very 
large values of Uoo are, however, permitted alternatives to standard general relativity. 

The decay of G in similar scalar-tensor theories would occur over time scales that are large 
compared to the Hubble time, and at late finite times the spacctimcs they predict would be 
observationally similar to de Sitter, although both G and H would be decaying with time. 
This is an interesting scenario. If space-time is asymptotically de Sitter, then H ^ Hq at 
late times. Since the expansion is accelerating, there exists an event horizon, and as is the case 
with a black hole horizon, the de Sitter horizon has a temperature: T = Hq/Itt = const 7^ 0. 
The presence of a non-zero minimum temperature in space-times that asymptote to de Sitter 
is important, as it means that even very rare thermal fluctuations are eventually expected to 
occur. For instance, the probability that a certain thermal fluctuation occurs goes as e~^l'^ 
for some action S which does not depend on T . The 4-volume of space-time grows as 
which in de Sitter spacetime grows as e^^°*. Thus, after a certain time i, the number of times 
that a particular thermal fluctuation is expected to have occurred is n ~ g-S/T-t-sm Provided 
T — ^ const, n will eventually be greater than unity, and n — ^ 00 asymptotically. A particularly 
extreme and topical example of this is the spontaneous emergence of "Boltzmann brains" in 
general -relativistic space-times that asymptote to de Sitter. Self-aware "observers" could 
emerge from the vacuum as a result of thermal fluctuations [20]. In an eternal, asymptoti- 
cally de Sitter, universe such "observers" shoTild vastly outnumber "ordinary" carbon-based 
observers such as ourselves. This conclusion is changed in scalar-tensor cosmologies. 

In scalar-tensor gravity, we have seen that theories in which either a; Wqo > 40, 000, or 
a; — ^ 00 as ^ — ^ 00, admit asymptotic solutions that would be observationally indistinguishable 
from de Sitter today, but in which G, — > at late times. Since G and H would, in such 
theories, decay only very slowly over a Hubble time, we still have T ~ H/Itt. The 3- volume 
of a spatial slice of constant t goes like a^/H^ oc g^v-^^^^ ^ and so the number of thermal 
fluctuations that are expected goes like: 

n ~ I* e-5/TW+3.-31nF^^ ^ J' j^^_^+3.-41nF^j ^ (^^-|f^+3.^ ^^^^ 

where, as above, F = H/Hq and y = Ina. 
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We have found that in ah theories where a; ^ oo or w ^ cooo > 40, 000 as — > cxd at late 
times, F ^ slower than e~'^y , where q < 5x 10~^. Even with this strong restriction, however, 
there is still a class of theories for which F — > faster than l/y as y ^ oo. Specifically, from 
Eqs. (33)-(34) we find that this will be the case if, as ^ ^ oo, grows more slowly than 
^1/2. Brans-Dicke theories with uj = const certainly satisfy this constraint. In these theories, 
as y ^ oo, the integrand in Eq. (45) would asymptotically tend to zero faster than e~^y for 
any r. Thus, n — > Hoo = const < oo at late times. If S is small enough, then no© ^ 1- In 
these cases, even though the universe would still be eternal, very rare thermal fluctuations 
like "Boltzmann brains" would not be expected to occur even once. 

This is another example (see also [21]), of how, if some or all of the traditional constants 
of Nature vary slowly with time, we cannot use current observations of the universe to make 
definitive statements about the expected behaviour of the universe in the far future. 

5 Conclusions 

We have investigated two features of the general behaviour of scalar-tensor gravity theories. 
Motivated by the need to understand the possible origins of de Sitter expansion in the early 
and late periods of the universe's history we have investigated how it can arise in general 
scalar-tensor gravity theories. We considered Friedmann universes filled with a mixture of a 
vacuum stress with equation of state p = —p = —po and other fluids which have total energy 
density pi = poK{y) and pressure pi, where y = Ina and limy^ooK(y) = liniy^oo Ky{y) = 
0,so the vacuum stress dominates at late times. In scalar-tensor theories with a coupling 
function uj{(p), we find that there is asymptotic approach to de Sitter expansion at late times, 
where cp — > (pooi provided the coupling function a;(^) diverges faster than \4>oo — ^l"^"*"*^ for all 
e > as — > ^ 0. This means that, for a given theory, specified by w(0), there must exist 
some ^oo e (0, oo) such that a; ^ oo and 0;'/^^+^ as (f) — > in order for cosmological 
solutions of the theory to approach de Sitter expansion at late times. This differs from the 
conditions required to establish a general relativity limit that has vanishing corrections to the 
weak-field PPN corrections to general relativity in the solar system: a; — *■ oo and co'/u^ 0. 

Brans-Dicke theory (a; = const) does not emit a de Sitter limit in the presence of p = — p 
stress. There is instead power-law inflation and G oc t^^. With this behaviour in mind 
we analysed the possible late time evolution of G{t) in the Friedmann cosmological models 
of scalar-tensor theories defined by an arbitrary Lo{(f)) and divided them into three classes 
depending upon the rate of decay of G with the expansion scale factor. The scenarios in 
which G decays over quickly (i.e. over a Hubble time or faster) would be difficult to realize in 
a manner that was compatible with solar system tests of gravity i.e. lo > 40, 000 today. The 
subset of theories with a physically viable slow decay were found to be particularly interesting 
because if, as is the case in Brans-Dicke theory, u ~ o{4>^/'^) as ^ — oo, then the expected 
number of extremely rare thermal fiuctuations, n(t), that occur within the visible universe, 
after a time t, would asymptote to a constant value. In general relativistic dc Sitter space- 
time, n{t) oo, which has led some to postulate that a typical 'observer' of our Universe 
would most likely have arisen out of the vacuum as a thermal fluctuation. The intrinsic 
probability of such a "Boltzmann brain" fluctuation is tiny, but since n{t) ^ oo in an eternal 
de Sitter universe, the number of such observers would grow without bound. In scalar-tensor 
theories, like Brans-Dicke, we have shown that the slow decay of G{t) can prevent this strange 
situation from occurring with any signiflcant probability. 
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